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ABSTRACT 



K. TUTTLE^, J. T. GREEN^ 



We present results from multi-epoch neutral hydrogen (HI) absorption observations of six bright pulsars 
with the Arecibo telescope. Moving through the interstellar medium (ISM) with transverse velocities 
of 10-150 AU yr~^, these pulsars have swept across 1-200 AU over the course of our experiment, 
allowing us to probe the existence and properties of the tiny scale atomic structure (TSAS) in the 
cold neutral medium (CNM). While most of the observed pulsars show no significant change in their HI 
absorption spectra, we have identified at least two clear TSAS-induced opacity variations in the direction 
of B1929-I-10. These observations require strong spatial inhomogeneities in either the TSAS clouds' 
physical properties themselves or else in the clouds' galactic distribution. While TSAS is occasionally 
detected on spatial scales down to 10 AU, it is too rare to be characterized by a spectrum of turbulent 
CNM fluctuations on scales of 10^"'^ AU, as previously suggested by some work. In the direction of 
B1929-l~10, an apparent correlation between TSAS and interstellar clouds inside the warm Local Bubble 
(LB) indicates that TSAS may be tracing the fragmentation of the LB wall via hydrodynamic instabilities. 
While similar fragmentation events occur frequently throughout the ISM, the warm medium surrounding 
these cold cloudlets induces a natural selection effect wherein small TSAS clouds evaporate quickly and 
are rare, while large clouds survive longer and become a general property of the ISM. 

Subject headings: ISM: clouds — ISM: structure — line; profiles 



1. INTRODUCTION 

For many years, observations of the diffuse interstellar 
medium (ISM) have traced a whole hierarchy of st r uctur es 
on spatial scales ^ 1 pc (cf. iDickev fc LockmanI (|1990[ )'). 
often attributed to the ISM turbulence. However, the ex- 
tremely small-scale end of this spectrum, reaching down 
to scales of tens to hundreds of AUs, has been largely un- 
explored. As the kinetic energy cascades from larger to 
smaller scales, it is expected that the turbulent spectrum 
will reach its end at so-called "dissipation scales." Under- 
standing of the dissipative properties of turbulence, such as 
its scale, efficiency, dissipation rate, and ubiquity, are im- 
portant since the delicate balance between energy injection 
and dissipation has a profound impact on the stellar ~ ISM 
recycling chain. However, many questions still await both 
theoretical understanding and observational constraints. 
For example, how is the balance between turbulent energy 
injection and dissipation maintained? What is the effi- 
ciency of the dissipation process? Does it vary with time 
and galactocentric radius? Is dissipation located only in 
pockets of "active turbulence", or it is ubiquitous? How 
do dissipative prope rties vary with gala ctic environments? 

Starting in 1976 (jDieter et al.l 119761 ) . structure in the 
diffuse ISM has been observed sporadically on scales down 
to 10^~^ AU. These findings caused much controversy as 
the observationally inferred properties were not in accord 
with the standard ISM picture. With a typical, observa- 
tionally inferred, HI volume density of ^ 10^ cm^^, and 
the thermal pressure of ^ 10^ cm~^ K (assuming tempera- 
ture of ~40-70 K), the AU-scale "tiny" features appeared 
significantly over-dense and over-pressured compared with 
the traditional cold neutral medium (CNM) clouds. 

Since the initial observations, the AU-sized struc- 



tures came to be rather frequently detected, with claims 
that up to 15% of the CNM could be in this form 
(JFrail et al.lll994[ ). There has been considerable progress 
in recent years in understanding of the ISM turbu- 
lence and the small-scale structure in the ISM, from 
both theo retical and observa t ional p erspectives [for re- 
views see lElmegreen fc Scala (|2004l ): IScalo fc Elmegreenl 
(2_0QJ)]- Nevertheless, it remains unclear whether the ob- 
served AU-scale features represent the dissipation scales 
of the ^ 1 pc spatial scale turbulence. At the same 
time, the observational picture has become even more 
complicated with improved telescope sensitivity, result- 
ing in the detection of new po pulations of possibly re- 
lated sub-parsec ISM clouds (IBraun fc Kanekad I2005I : 
IStanimirovic fc Heilesll2005t iBegum et al.ll2Q10[ ). The "low 
column density clouds" observed by Stanimirovic fc Heiles 
(2005) have a peak optical depth of only 10~^ to 10"^ and 
an inferred size of 500-5000 AU. 

Motivated by the recent theoretical and observational 
efforts in understanding the nature and origin of the tiny- 
scale ISM structure, we have undertaken multi-epoch ob- 
servations of the HI absorption against a set of bright pul- 
sars. Because of pulsars' relatively high proper motion 
and transverse speeds (typically ■^ 10 — 10^ AU yr~-^), 
pulsar HI absorption profiles obtained at different epochs 
sample CNM structure on AU spatial scales. In addition, 
the pulsed nature of pulsars' emission allows spectra to be 
obtained on and off source without moving the telescope, 
therefore sampling both emission and absorption along al- 
most exactly the same line of sight. IStanimirovic et al.l 
(2003) summarized preliminary results from this project 
for three pulsars, B0823-H26, B1133+16 and B2016+28. 
Here we present results for all pulsars with full analysis. 
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The structure of this paper is organized in the following 
way. Section [2] provides a short summary of the main ob- 
servational and theoretical results of the AU-scale struc- 
ture in neutral gas. In Section [3] we summarize our ob- 
serving and data processing strategies. Results from our 
investigation of changes in absorption profiles' integrated 
optical depth properties (c.f. equivalent widths) are pre- 
sented in Section |4l while in Section [5] we search for HI 
absorption feature optical depth variability. To facilitate 
further data analysis we derive the spin temperature and 
HI column density in the direction of all pulsars in Sec- 
tion [Bl In Section 7, we investigate in detail the properties 
of the ISM along lines of sight toward our pulsars. In 
Section |8] we study the properties of the most significant 
AU-scale features found in this study. We discuss our re- 
sults and compare them with several theoretical models in 
Section [HI and summarize our findings in Section [TU] 

2. BACKGROUND 

While structure on AU scales has been observed in 
neutral, ionized, and molecular flavors of diffuse gas, it 
is still not clear whether all these features trace a sin- 
gle or multiple phenomena. Our work focuses only on 
the tiny-scale atomic structure (TSAS; Heiles 1997) and 
we summarize results from two different observational 
approaches used to study TSAS in the CNM below in 
i j2.1l and ^2.21 as well as several theoretical considera- 
tions of the TSAS phenomenon in iJ2.3l In addition 
to observational methods discussed here, TSAS in the 
CNM has been also observed through temporal (prob- 
ing scales of ^ 1 to tens of AUs) and spatial (probing 
scales of a few thousands of A U) variability of optical ab- 
sorption lines (Mever fc Blades 1996; Andrews et al. 20Ql|; 
ICrawford et al.l l2000t iLauroesch fc Meveij I2003D . The 
reader is refer red to excellent articles in the ASP Conf. Se- 
ries Vol. 365 ([Haverkorn fc G oss 2007) for complementary 
details on optical absorption line studies, or the Galactic 
tiny-scale structure in the ionized and molecular gas. Fi- 
nally, there are also some indications that even damped 
Lyman a absorbing systems at high redshift may con- 
tain small-scale HI structure on scales of ~ 10 — 100 AU 
(Kanekar fc Chengalur 2001). 

2.1. Observations of TSAS via spatial mapping of HI 
absorption line profiles across extended background 
sources 

Single baseline very long baseline interferometric 
VLB I) observations of the quasar 3C147 bv lDieter et al.l 
1976[ ) were the first to infer the existence of a tiny-scale 
atomic cloud, having a size of 70 AU and HI volume den- 
sity of 10^ cm~^, based on variations in absorption line 
profiles in the directions of separa te source comp o nents . 
Subsequent VLBI observations by [Diamond et al.l (1989) 
found evidence for 25-AU HI clouds in directions toward 
3C138, 3C380, and 3C147. Particularly large variations 
were noticed in the case of 3C138. 

The first images of HI optical depth distribution 
in the directio n of ex tragalactic sources were obtained 
by iDavis et al.l (| 19961) toward 3C138 and 3C147, using 
the MERLIN arra y and the European VLBI Network. 
iFaison et all (| 19981 ) and iFaison fc GossI (|2001l ) used the 
Very Long Baseline Array (VLBA) to image HI absorption 



toward seven sources. Significant variations in HI optical 
depth with At of 0.1 and ^ 1, respectively, were found 
toward only two sources, 3C147 a nd 3C138. 

The latest results on 3C138 bv lBrogan et al.l (|2005l ). in 
a three-epoch series of observations, show clear evidence 
for spatial variations of At ~ 0.5 on scales of 25 AU. In 
addition, temporal changes in HI optical depth have been 
found over a period of 7 years, wit h implied tra n sverse 
velocities of order 20 km s~^ (see iBrogan et al.l (|2007f) 
for a thorough review of these measurements). Some 10% 
of pixels in their optical depth images have At > 5-cr. 
Assuming that most of the CNM along the line of sight 
contains small-scale structure, they estimated the volume 
filling factor of about 1%. They proposed that TSAS is 
ubiquitous, and that the lower measured levels of varia- 
tions in the other cases result from a selection effect due 
to incomplete sampling of the relevant angular and spatial 
scales. 

Recentlv. lGoss et al.l ()2008[ ) used MERLIN observations 
of 3C161 and found At ~1.0-1.5 on scales of ~ 500 AU. 
A slightly smaller level of variations (typically At ■^ 0.2) 
was found in the c ase of 3C147, on scales of about 10 AU 
(jLazio et al.ll2009[ ). 

2.2. Observations of TSAS via time variability of HI 
absorption profiles against pulsars 

In the late 1980s, sufficiently accurate and repeated 
measurements of HI absorption profiles against pulsars be- 
gan to be made, and it was noticed that pulsar ISM spec- 
tra changed over time in some cases, suggesting inhomo- 
.eneit ies in the intervening gas. For example, Clifton et al.l 
1988() found that the HI absorption spectrum of PSR 
B1821-h05 changed significantly between 1981 and 1988, 
with the appearance at the latter epoch of a previously 
unobserv ed featu re with t ^ 2 and Aw r^ 1 km s~^. 
JDeshpande et al.l ([1992) showed that between ~ 1976 and 
1981, HI absorption toward B1154-62 did not change sig- 
nificantly; while toward B1557-50, a variation with At ^ 1 
was interpreted as a cloud of size in th e 1000 AU range. 

The early pulsar HI results inspired iFrail et aLl (|1994[) 
to undertake a dedicated multi-epoch pulsar HI absorp- 
tion experiment at Arecibo. Six pulsars were observed 
at three epochs, with time baselines ranging from 0.7-1.7 
yr. These authors reported the presence of pervasive vari- 
ations with At r^ 0.03 — 0.7, and associated HI column 
densities of 10^^ - 5 x 10^° cm~^; on scales of 5-100 AU. 
They indicated that 10-15% of cold HI is in the tiny struc- 
tures, and additionally detected a correlation between ab- 
sorption equivalent width variations and equivalent width 
itself. These results appeared to buttress the earlier VLBI 
findings of TSAS, and they provided a strong impetus for 
further experimental and theoretical work. 

The recent era of pulsar TSAS experimen ts began with 
the Parkes observations of iJohnston et al.l (2003' ) . Sur- 
prisingly, these investigators found no significant optical 
depth variations in their three-epoch, 2.5-yr observations 
of three pulsars, and were able to place an upper limit on 
column density variations t oward PSR B1641-45 of 10^^ 
cm^^, significantly below the lFrail et al.l (|1994[ ) detections. 
They showed that the earlier experiment did not fully ac- 
count for the large increases in noise in absorption spec- 
tra at the line frequency, so that some of the apparently 
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significant variations actually were not. Whi l e no s ignif- 
icant variations were seen by iJohnston et alj (l2003l ) dur- 
ing the 2.5 years of their TSAS experiment, they detected 
variations in the spectrum of PSR B1557-50 when com- 
pared with measurements made five years earlier. This 
is the same pulsar who se spectrum was noted earlier by 
iDeshpande et al.l ()1992[ ) to vary on similar timescales in 
the late 1970s. In combining the results from fou r mea - 
surements over twenty-five years, IJohnston et al.l (|2003f ) 
concluded that the cloud causing the variations is ~ 1000 
A U in size, with a de nsity of ^ 10^ cm~^. 

iMinter et al.l (|2005[ ) performed a particularly exhaustive 
TSAS study on PSR B0329-I-54, which is very bright and 
almost circumpolar at the Green Bank Telescope. The in- 
vestigators observed the pulsar continuously for up to 20 
hr in eighteen observing sessions over a period of 1.3 yr. 
They detected no HI optical depth variations (At < 0.026 
in most cases) for pulsar transverse offsets ranging from 
0.005-25 AU. 

We and our colleagues published early results 
from a new set o f multi-epoch Arecibo observations 
(jStanimirovic et aLl 2003). Based on non-detections of 
TSAS in the direction of three pulsars (B0823+26, 
B1133-hl6, and B2016-h28), we suggested that TSAS is 
not ubiquitous in the ISM and could be related to isolated 
events. The current work elaborates these themes. 

2.3. Theoretical considerations of TSAS 

Several explanations were prop osed to reconc ile the ob- 
servations of TSAS with theory. iHeilesI ()1997( ) proposed 
that TSAS features could be curved filaments and/or 
sheets that happen to be aligned along our line-of-sight. 
By combining low temperature and a geometrical elonga- 
tion, the apparently problematic volume density and ther- 
mal pressure can be significantly reduced. We discuss this 
m odel in som e detai l in Section IHTT] 

iDeshpandg (|2000D suggested that TSAS represents the 
tail-end of the turbulent ISM spectrum which exists on 
larger spatial scales, and showed that the observed opti- 
cal depth differences are consistent with a single power- 
law description of the HI optical depth distribution as a 
function of spatial scale. It was also pointed out that the 
over-density of TSAS is due to a misinterpretation of ob- 
servations: the observed variations in optical depth sample 
the square-root of the structure function of the HI opti- 
cal depth, not directly the power spectrum of optical depth 
fluctuations. If the power spectrum of the HI optical depth 
distribution in the dire ction of Gas A, with a measured 
slope of 2.75 (Deshpan de et al.ll2000[ ) over a range of spa- 
tial scales from 0.02 to 4 pc, is extrapolated to AU scales 
then At = 0.2-0.4 is expected on scales of 50-100 AU. We 
note that this slope was confirmed re cently in an indepen- 
d ent experimen t bv lRov et al.l ()2010D . 

iGwinnI (|2001h proposed that optical depth fluctuations 
seen in multi-epoch pulsar observations are a scintillation 
phenomenon combined with a vel ocity gradient acro ss the 
absorbing HI. While in the case of lDeshpandd (J2000() opti- 
cal depth variations are expected to increase with the size 
of structure, the scintillation phenomenon predicts maxi- 
mum variations on the very small spatial scales probed by 

•^ The Arecibo Observatory is part of the National Astronomy and Ionosphere Center, operated by Cornell University under a cooperative 
agreement with the National Science Foundation. 



interstellar scintillatio n. However, dedica ted observations 
of PSR J1456-6843 bv lGwinn etall (I2007D failed to detect 
any scintillation effects on the HI absorption line. 

One of the main problems regarding TSAS is its over- 
pressure relative to the traditional thermal pressure of 
the ISM, which suggests that TSAS can not persist 
nor be pervasive. However, recent numerical simula- 
tions of the dynamic and turbulent ISM find substantial 
departures from pressure equilibrium. Instead of well- 
defined thermal phases, gas with a wide r ange of densities 
and temperatures is found in sim ulations (|Mac Low et al.l 
l2005t iHennebeUe fc AuditI [20071 ). Such results are sup- 
ported by observations of the GI fine-structure lines by 
iJenkins fc Trippl (|200l[ ) who found evidence for the exis- 
tence of over-pressured gas {P/k ^ 10^ cm~'^ K). While 
these studies suggest that high-pressure TSAS could be a 
natural product of the highly dynamic ISM, the fraction of 
high-pressure gas in the ISM, as well as the physical pro- 
cesses responsible for their formation, vary hugely across 
simulations. The only way to constrain TSAS properties 
and abundance is through sensitive observations. 

3. OBSERVATIONS AND DATA PROCESSING 

3.1. The observing and basic data processing procedure 

We used the Arecibo telescopqfl and the Ga l tech B ase- 
band Recorder (GBR) backend (jJenet et al.l I1997D . to 
measure HI absorption spectra of six pulsars at multiple 
epochs. We chose to observe the same sources as Frail et 
al. (1994) in order to enhance the number of available 
time baselines for comparison. Table [T] lists several basic 
parameters for our targets: Galactic coordinates (l/b), the 
dispersion measure (DM) , continuum equivalent flux den- 
sity at 20 cm (S20), distance (D), the transverse velocity 
(Vt), the transverse distance traveled in a year (L_l), and 
the pulsar-on antenna temperature (Tpsr). The tabulated 
flux density and TpsR represent average values obtained 
from long integrations, whereas the actual pulsar flux den- 
sity and TpsR vary on timescales from several minutes to 
several hours due to interstellar scintillation. Distances 
for all six pulsars (column 5) were determined from inter- 
ferometric parallax measurements. The pulsar-on antenna 
temperature Tpsr = GS20P/W50, with the L-band gain 
G = 7.5 X 10-3 K mJy-i (Heiles et al. 2001), and pulsar 
period P and 50% pulse width W50 from the ATNF pulsar 
Gatalogue (Manchester et al. 2005). 

The observing and basic data proces s ing p roce- 
dures were des c ribed in IStanimirovic et ahl (J2003D and 
IWeisberg et al.l (|2008[ ). Rela tive to t he previous large 
multi-epoch experiment by F rail et al.l (1994), our four- 
level spectrometer is less prone to systematic errors when 
observing these strong and highly- variable sources. We 
had four observing sessions, referred to below as SI, S2, 
S3, S4; for epochs 2000.62, 2000.95, 2001.70, and 2001.86, 
respectively. The typical integration time per pulsar was 
about 8 hours for each of four sessions. 

The raw complex voltage samples recorded with the 
GBR were processed at Galtech's Genter for Advanced 
Gomputation and Research to obtain a data cube of pulsar 
intensity as a function of pulsar rotational phase and radio 
frequency. For each scan, portions of the data cube during 
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Fig. 1.— Emission (top) and absorption (bottom) profiles in tlie direction of B0823+26, B1133+16, B1737+13, and B2016+28 obtained at 
tlie first opocii (August 2000). The emission spectra were scaled to match the brightness temperature from the Leiden-Dwingeloo survey. 
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Fig. 2. — Emission and absorption profiles in the direction of B0540+23 and B1929+10, obtained at the first epoch (August 2000). The 
middle panels show the absorption profiles before baseline fitting for B0540+23, and before ghost removal (see i|3.2l l for B1929+10. In the 
case of B0540+23, a large ripple is obvious in the middle panel, caused by the interstellar scintillation. For B1929+10, a faint inverted copy 
of the H I emission spectrum is noticeable especially around 20 km s~^. The bottom panels show the absorption spectra after baseline and 
ghost correction, respectively. 



Table 1 

Basic parameters of pulsars observed in this study, tabulated or derive d using the reference given in the last 

column and the atnf pulsar catalogue ( manchester et al.|[20b5i l. 



PSR 


l/b 


DM 


(B1950) 


deg 


(pccin" 


B0540+23 


184.36/-3.32 


77.7 


B0823+26 


196.96/31.74 


19.5 


B1133+16 


241.90/69.19 


4.8 


B1737+13 


37.08/21.68 


48.9 


B1929+10 


47.38/-3.88 


3.2 


B2016+28 


68.10/-3.98 


14.2 



S20 
(mJy) 



D 

(kpc) 



Vt 
(km s" 



') 



(AU yr- 



TpSR 

(K) 



Reference 



9 


3.5 


377 


10 


0.4 


192 


32 


0.4 


631 


3.9 


4.7 


672 


36 


0.3 


177 


30 


1.0 


30 



80 
39 
135 
142 
37 
7 



2.7 
6.9 
9.0 
0.7 
8.3 
8.4 



Harrison et al. fl993) 
Gwinn et al. (1986) 
Briskcn et al. C2002) 
Bris ken et al. (2003) 
Chatteriee et al. (2004) 
Brisken et al. (2002) 



Pulsar Measurements of Tiny-Scale Atomic Structure 



and off the pulsar pulse were accumulated to obtain the 
'pulsar-on' and 'pulsar-off' spectra. The pulsar absorption 
spectrum represents the 'pulsar-on' - 'pulsar-off' spectrum 
for each scan, whose baseline was then flattened by per- 
forming frequency switching and normalizing by the off- 
line pulsar intensity. In accumulating pulsar absorption 
spectra from different ^ 2— min scans, each spectrum was 
given a weight proportional to TpsR, where TpgR is the 
antenna temperature of the pulsar. 

The frequency-switched 'pulsar-off' spectra correspond 
to the H I emission spectra in pulsar directions. As no 
temperature calibration was performed during the obser- 
vations, we have scaled these spectra to match the H i 
brightness temperature from the Leiden-Dwingeloo survey 
(jBurton fc Hartmannlll994 ). Final absorption and emis- 
sion spectra for all pulsars, from the first observing epoch 
(2000.62), are shown in Figures [1] and [2J Our spectrome- 
ter, consisting of 4096 frequency channels across 10 MHz 
total bandwidth, has a velocity channel spacing of 0.52- 



km 



However, adjacent spectral channels are not in- 



dependent, which can lead to distortio n of narrowband 
signa ls such as the ones we are studying (|Rohlfs fc Wilson! 
l2004f ). To ameliorate this problem, we Hanning smoothed 
all processed absorption spectra, yielding a final velocity 
resolution of 1.04 km s~^. Note that Hanning smoothing 
was not applied in our preliminary analyses [Stanimirovic 
et al. (2003); Weisberg & Stanimirovic (2007)]. 

3.2. Baseline fitting of absorption spectra 

To flatten the baseline of the final absorption spectra 
we fitted a polynomial function to the unabsorbed por- 
tions of the spectrum. For pulsars B0823+26, B1133+16, 
B1737-hl3 and B2016-H28 generally a 1st or 2nd order 
polynomial function was fitted and then the full spectrum 
was divided by the fitted function. A higher order polyno- 
mial was required for the 3rd session data for B2016+28 
due to ripples in the baseline structure, most likely cased 
by interstellar scintillation. 

PSR B0540+23 data were heavily affected by interstellar 
scintillation. T his pulsar's sc intillation bandwidth at 1.4 
GHz is 0.3 MHz (|Gwinnll200ll) . which resulted in numerous 
ripples across our band (see Figure [U left middle panel). 
We needed to fit and divide by a 15th order polynomial 
to fiatten the baseline for this pulsar (see bottom panel in 
Figure m. However, as emphasized in Section 4.1, we do 
not use this pulsar's data for our analysis due to large scin- 
tillation effects which presumably led to the unphysically 
(///o < 0) deep absorption features. Noise and calibration 
errors in the pulsar absorption spectrum are magnified by 
a factor of e"^ when the absorption spectrum is expressed in 
terms of the ultimately desired optical depth r. With this 
pulsar's optical depth near infinity, these errors become so 
large as to render our measurements unusable. 

B1929-I-10 data, as shown in Figure p (right, mid- 
dle p anel), suffer from the 'ghost effect' (jWeisberg et al.l 
Il980| ). The main absorption feature is at 5 km s~^. How- 
ever several additional features are visible from 20 to 80 
km s~^ which represent a false image of the pulsar-on 
emission spectrum. The usual reason for the ghosts is an 
incorrect scaling of total power received by the spectrome- 
ter in the case of strong and variable sources. As the ghost 
is a faint copy of the pulsar-on emission spectrum in the 



pulsar absorption spectrum, we can exorcize the ghost by 
fitting for and removing the former from the latter. 

To do so we have performed a least-squares fit of the 
I(v)/Io absorption spectrum to a ghost of fitted ampli- 
tude A, and a fitted linear baseline slope B plus constant 
C, all to spectral channels lacking true absorption: 

I{v)/Io = ATpsn-on{v) + Bv + C, (1) 

where Tpsn-oniv) is the measured pulsar-on spectrum. 
The exorcized absorption spectrum is then: 

(/(^)//o)corrected = [I{v)/Io " ^rpSR-on(^) " Bv]/C. (2) 

The exorcized PSR B1929-I-10 absorption spectrum is 
shown in Figure [21 right bottom panel. The fitted param- 
eters are listed in Table [2] 

3.3. Noise envelopes 

For each absorption spectrum we have calculated an ex- 
pected noise envelope which will be used to assess the 
significance of any observed temporal variations. The ex- 
pected noise level AT{v) is proportional to the overall sys- 
tem temperature Tsys{v), which has the following contri- 
butions. (1) The receiver temperature Trx ^ 25-30 K. (2) 
The sky background, which is typically of o rder of 1-2 K, 
based on the all-sky survey at 408 MHz by iHaslam et al.l 
(|1982[ ) and after scaling to 1.4 GHz using a spectral in- 
dex of —2.6. (3) The contribution from the HI emission 
which can be very significant and which varies with veloc- 
ity. For example the rms noise on-line is three times higher 
than the rms noise off-line for the case of B0540-I-23. (4) 
At Arecibo, the pulsar continuum emission itself is also 
sometimes sufficiently strong to provide the principal con- 
tribution to the system temperature. 

Since an observed pulsar-on spectrum is identically the 
sum of these four contributions, the expected noise in the 
pulsar-on spectrum should be proportional to the pulsar- 
on spectrum. Furthermore, owing to pulsar-off integra- 
tion times being '^ 10 — 100 times longer than pulsar-on 
integration times, the pulsar-off spectrum adds negligible 
additional noise to the absorption spectrum. Hence our 
expected noise envelopes (shown in Figures [S HlOp are gen- 
erated by scaling the observed pulsar-on spectra by a mul- 
tiplicative factor that matches them to the observed noise 
in the off-line portions of the absorption spectra. We have 
then propagated uncertainties associated with baseline fit- 
ting and ghost removal procedures to derive the final ex- 
pected noise envelopes. To estimate significance envelopes 
for difference spectra shown in Section [5] (Figures 3-8) or 
for the equivalent width analysis in Section 21 appropri- 
ate noise envelopes for individual absorption profiles were 
added in quadrature. 

To investigate the properties of noise in emission-free 
channels, we have plotted the histogram of data points lo- 
cated off-line for each difference spectrum (not shown in 
the paper) . Most histograms are consistent with Gaussian 
statistics with a standard deviation in agreement or close 
to our derived 1-a noise level. In several cases there is a 
small discrepancy between the measured off-line noise and 
derived 1-a noise envelopes. This is not surprising after 
the propagation of errors due to baseline, and especially 
ghost, fitting. For at least three of our pulsars we have 
achieved excellent sensitivity. We discuss the rms on-line 
noise level of the final absorption spectra in Section [51 
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Table 2 
Fitted parameters for the 'ghost' excision in the case of B1929+10 (Section 3.2). 



Session 



A 



(10- 



B 



C 



1 


-4.3 ±0.4 


-3.0 ±0.6 


1.148 ±0.003 


2 


-9.6 ±0.4 


6.1 ±0.3 


1.147 ±0.002 


3 


-7.5 ±0.4 


4.8 ±0.3 


1.128 ±0.002 


4 


-14.9 ±0.7 


-0.3 ±0.3 


1.144 ±0.002 



4. EQUIVALENT WIDTH VARIATIONS 

Our main goal is to search for variations in the HI ab- 
sorption profiles obtained at different observing epochs. 
We have employed two methods: (i) in this section we look 
for changes in integrated properties such as the equivalent 
width, and (ii) in Section [S] we investigate spectra differ- 
enced between each pair of observing epochs. 

We have measured the equivalent width 



EW 



rdv 



HI/Io)dv, 



(3) 



in units of km s~^, for each observing epoch. The opti- 
cal depth profiles were integrated over the 'on-line' veloc- 
ity regions. To avoid low signal/noise parts of absorption 
spectra, as well as line wings which may be affected by 
small baseline imperfections, we have defined the 'on-line' 
regions as roughly encompassing velocities with t{v) >2-ct. 
For each pulsar a single region was used for all observing 
epochs: [-10, 15] km s'^ for B0540+23, [2,7] km s"! for 
B0823+26, [-7,0] km s'^ for B1133+16, [2,6] km s^^ for 
B1737+13, [-3,12] km s^^ for B1929+10 and [-2,21] 
km s^^ for B2016+28. In addition, by averaging optical 
depth profiles for all epochs we have estimated the mean 
equivalent width for each pulsar (EW). 

To compare EW from one epoch to the next we have cal- 
culated AEW(Ai) = [EW(ii)-EW(i2)]- These values, as 
a function of time baselines At (bottom axis) or traversed 
spatial scales I (shown on the top axis), are shown for all 
pulsars in Figure [3l B0540+23 was excluded from this 
analysis as its optical depth profile is uncertain around 
km s~^. To convert from time intervals At to spatial scales 
I, we use / = At X L±, where L± is provided in Table 1. 
Our spatial scales given below should be multiplied by an 
(unknown) fractional distance factor / = -DhiZ-Dpsr to 
account for the fact that absorbing HI is located between 
us and the pulsa r. Note that some au thors have arbitrarily 
fixed /; e.g. for iMinter et"aLl (|2005[ ) / = 0.5. (When we 
assemble various workers' results in this paper, we redefine 
all such values onto our scale.) 

As shown in Figure [3l in the case of B0823+26, 
B1133+16, B1737+13 and B2016+28 all variations in 
AEW are essentially within 2-cr and are most likely not 
significant. However, in the case of B1929+10, baselines 
of 0.33 and 1.24 yrs show changes in AEW at a 3.3-cr and 
3.5-cr level, respectively. In addition, the 0.16 yr baseline 
shows a variation at a 2.8-ct level, while the 0.75 yr base- 
line shows a 2.5-cr change. To quantify the significance of 
AEW variability we test the null hypothesis of AEW being 
constant with time. We fit a constant AEW model to mea- 
surements in Figure [3] and find xV^of =0.4 (B0823+26), 



0.8 (B1133+16), 1.3 (B1737+13), 6.0 (B1929+10), and 
1.2 (B2016+28). We conclude that AEW changes for 
B1929+10 can not be fit with a constant model, while 
for all other pulsars no significant variability of AEW is 
observed. This AEW variability, as further analyzed in 
Section 7.1, can be interpreted as being due to TSAS. 

In Figure |4] we plot AEW, separately for individual 
pulsars, but as a function of EW. This plot can be di- 
rectly compared with Figure 4 from Frail et al. (1994). 
Our results for B0823+26, B1133+16 and B1929+10 agree 
roughly with Frail et al.'s, | AEW [< 0.1. However AEW 
for B1737+13 and B2016+28 found in Frail et al. (1994) 
were significantly larger than what we find. A particu- 
larly striking difference is the case of B2016+10, where 
Frail et al. (1994) found AEW - 1-5, while we find only 
AEW ~ 0.1-0.3. We believe that at least some of the 
variations Frail et al. saw were due to small calibration 
errors. Looking at our whole sample, the fractional vari- 
ation AEW/EW spans a range from 0% to 11%, with a 
median of 4%. This is significantly smaller from a median 
value of 13% found in Frail et al. (1994). 

As a conclusion, we find significant variations in AEW 
only for B1929+10, with two largest changes being at a 
fractional level of 11%. The frequency of variations and 
the level of variations (as well as upper limits) are sig- 
nificantly smaller than what was found previously for the 
same sources by Frail e t al. (1994 ) . Ou r results do not 
rule out the finding by iFrail et all ()1994 D that [ AEW [ 
increases with EW. However, if this relationship exists, it 
is at a much smaller level than suggested by Frail et al. 

5. variations from difference spectra 

Considering that TSAS may have temperature as low 
as 15 K (Heiles 1997), the change in absorption profiles 
may occur over <^ 1 km s^^ and appear in only one or 
two velocity channels in our spectra. Such narrow changes 
may get washed out when calculating the equivalent width. 
We therefore in this section compare all pairs of absorp- 
tion profiles for a given pulsar, and investigate absorption 
spectrum differences, A(///o) ~ At. With four observing 
epochs we have six difference spectra (Figures [5] to [TU|). 
Overlaid atop each difference spectrum is a ±2-cr signifi- 
cance envelope calculated from noise envelopes discussed 
in Section [531 

In summary, we see changes in difference spectra at a 
3-cr level only in the case of B1929+10 and only for the two 
baselines where we also find changes in AEW at a > 3-cr 
level. We note that B1133+16 has one feature above 3- 
(T (s4-s3 baseline) at a velocity of —30 km s~^, however 
this likely comes from RFI as is seen in emission in the 
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Fig. 3. — Temporal change in the equivalent width between two epochs, AEW(At) = [EW{t2) — EW(4i)]. Plotted error bars represent 
±2-(T uncertainties. The top x-axis on each plot shows the corresponding spatial scales, I = At X Lj^ . 
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Fig. 4. — The absolute value of AEW as a function of the mean EW for each pulsar, EW. Plotted error bars represent it 1-cr uncertainties. 
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absorption spectrum from session 3. Similarly, B1737+13 
has a significant variation at a velocity of ~ 35 km s^^on 3 
baselines, all involving the session 1 absorption spectrum; 
this corresponds to an (unphysical) emission feature and 
is most likely due to RFI. 

For three of our pulsars (B0823+23, B1133+16 and 
B1929+10) we have excellent sensitivity for TSAS exper- 
iments. At ^ 0.01-0.02 (2-CT in the on-line portion of 
the spectrum). Only one recent experiment (Minter at 
al. 2005) reached similar sensitivity but primarily over 
smaller spatial scales. For B1737-I-13 and B2016+28 we 
have reached 2-cr sensitivity of 0.1-0.2. Although this is 
slightly higher, several pulsar and VLBA TSAS detections 
have been claimed at this level. We provide below brief 
comparison between our observations for individual pul- 
sars and previous TSAS experiments. 

5.1. Difference spectra of individual sources 

5.1.1. PSR B0540+23 

The corresponding tranverse distances traveled by 
B0540-I-23 during the course of our observations are 13, 
26, 60, 73, 86 and 100 AUs. Frail et al. (1994) concluded 
that there were clear signs of significant optical depth vari- 
ations in some absorption components for this pulsar. In 
fact both our and Frail et al.'s spectra show (I/Iq) slightly 
less than zero near km s^^ as well as scintillation rip- 
ples. As already discussed in Section 3.2, we do not use 
data for this pulsar for our further analyses and do not 
show difference spectra. 

5.1.2. PSR B0823+26 (Figure\^ 

With L±_ = 39 AU yr^^, the time intervals covered by 
our experiment translate to spatial scales of 6, 13, 29, 35, 
42 and 48 AUs. There are no data points with Ar above 
3-cr in difference spectra for all six time baselines. We 
conclude that no significant change in absorption profiles 
is found down to a At level of 0.03 (a 2-cr level on H i line). 
Frail et al. (1994) detected only marginal changes over 0.6 
yr, but reported significant variations of At ~ 0.07 in 1.1 
yr. We do not confirm this reported variation on any of 
our temporal baselines although our 2-a sensitivity level 
on all baselines is at least two times smaller than the level 
of variations fo und by Frail et al. Q994). We discussed 
this further in IStanimirovic et aLl (|2003.') and concluded 
that calibration and noise issues affected the Frail et al. 
analyses. 

5.1.3. PSR B1133+16 (Figure\^ 

As this pulsar has L± — 135 AU yr~^ a wide range of 
spatial scales is probed by our observations, from 20 AU 
to 170 AU. No variability at ^ 3-cr was found and we con- 
clude that no significant variations were detected in the 
direction of this pulsar down to a At level of about 0.02. 
Frail et al. (1994) reported At ^ 0.04 on a 1.1 yr base- 
line. On all of our difference spectra (time baselines of 
0.33, 0.75 and 1.08 yr) our 2-cr sensitivity is half the level 
of variations found in Frail at al. (1994), At ^ 0.02 at 
most. We note that no variability was detected at this 
level in any of later variability experiments. 

5.1.4. PSR 31737+13 (Figure^ 



With ix — 142 AU yr~^ the range of probed spatial 
scales is from 20 AU to 175 AU. No significant varia- 
tions are found. Frail et al. (1994) reported At — 0.5 
at 4 km s^^(which corresponds to A{I/Io) ^ 0.1 in this 
relatively deep line) over 0.65 yr. Our 2-cr sensitivity is 
A(///o) ~ 0.2 for half of baselines, and 0.3 for the rest. 
Several TSAS features were detected at a similar level in 
recent experiments, both pulsar and interferometric. 

5.1.5. PSR B 1929+ 10 (Figure\^ 

With L± = 37 AU yr^^, the range of probed spatial 
scales is very narrow, from 6 to 45 AU. The absorption 
profile for this pulsar consists of at least three different 
velocity components (at —1.5, 5 and 8 km s"^). The com- 
ponent at 5 km s~^, which is the one with the deepest HI 
absorption, is the only one exhibiting significant changes. 
Difference plots show line variations at 3.2 and 3.0-cr lev- 
els (At =0.025) for time baselines of 0.33 and 1.24 years. 
We consider these detections to be significant, especially 
as the same time baselines showed significant variations in 
EW. This level of variability is the lowest ever detected in 
TSAS experiments. Both detected difference features are 
sharp and appear in the same single velocity channel at 5 
km s~^ over different epochs. We further measure and an- 
alyze properties of the intervening absorbing clouds traced 
by these detected features in Section[Sl Baselines of 0.16, 
0.75, 0.91, and 1.08 years have an upper-limit in At of 
0.012-0.018 at this velocity. Frail et al. noticed variations 
of up to At '^ 0.03 for this pulsar, with the largest fiuctu- 
ations being at a velocity of '^ 8 km s"^, where we observe 
no significant variations. 

As this is the only pulsar in our sample to exhibit signif- 
icant absorption variations, we have scrutinized the data 
acquisition and analysis procedure with extra care. For 
example, a scaling error of approximately 10% in pulsar 
strength in a single epoch could result in an erroneous 
rescaling of optical depths by a similar amount, which 
would lead to a spurious detection of variation when there 
is actually none. For example, in Stanimirovic et al (2003), 
we pointed out an apparent case of this type of problem 
in Frail et al's spectra of PSR B2016+28. However, it 
appears that the more advanced CBR spectrometer is in- 
deed able to measure spectra more accurately, as we an- 
ticipated. For example, we see no evidence of variations in 
our experiment with PSR B2016+28 nor with B1133+16, 
the other particularly bright pulsars in our sample. (Cal- 
ibration errors would be expected to be most serious in 
the brighter pulsars, which can easily increase the system 
temperature by factors of more than two.) 

It is notable that B1929-I-10 is the only pulsar in our 
sample from which we removed a ghost. Then an obvious 
question is whether our detected variability could be an 
artefact of the ghost fitting procedure. We have scruti- 
nized our reduction procedures and remain confident that 
the variability is real. Spectra for B1929-I-10 were reduced 
in the same way as for all other sources. Even the ghost 
fitting procedure was done in exactly the same way (using 
exactly the same IDL routines) as the polynomial base- 
line fitting for other pulsars. The derived error envelopes 
reflect additional uncertainties due to this fitting process. 

In Figure [S] we show all four absorption profiles for 
this pulsar before (top) and after (bottom) ghost removal. 
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Fig. 5. — Difference spectra for B0823-I-26 derived from four different observing epochs. Dashed Unes show the expected ±2- and 3-cr noise 
envelopes, whose shape reflects the HI emission in the beam. The vertical scale is in the original units of the absorption spectra, A(//7o), 
which, for small r, is also equal to A(r). The 'on-line' region is shown with a horizontal bar. 
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Fig. 6. — Difference spectra for B11334-16 derived from four different observing epochs. See Fig. [5]for further details. 
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Fig. 7. — Difference spectra for B1737+f3 from four different epochs. See Fig.[5]for furtficr details. 




Fig. 8. — Difference spectra for B1929+f0 derived from four different observing epoclis. Dashed lines show the expected ±2- and 3-cr noise 
envelopes, whose shape reflects the HI emission in the beam plus uncertainties associated with ghost excision (See t|3.2ll . 
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Fig. 9. — HI absorption spectra for B1929+10 from all four epochs before (top) and after (bottom) ghost removal. Dashed lines show 
the expected ±2- and 3-cr noise envelopes for the 1st epoch spectrum The HI emission spectrum, which was exorcized from the displayed 
absorption spectra (see text), is shown as a faint black line at normalized intensity ~ 0.92. 



Only the deepest velocity component at 5 km s~^ shows 
variability (both before and after ghost removal), while 
components at —1.5 and 8 km s~^ agree extremely well 
between epochs. If there were some imperfections related 
either to a scaling problem or to ghost fitting errors, they 
would affect all features. We have also examined the es- 
timated coefficients used in equation (2). In the case of 
the 1st session spectrum, for example, A = —4.3 ± 0.4. A 
smaller value for A would result in a larger correction and 
would bring the 1st session absorption spectrum closer to 
spectra from other sessions. However, even by taking val- 
ues within ±2 — 3a of the best-fit estimate, the observed 
variations in the difference spectra would still stay sig- 
nificant (above 3-ct). The reason for this is that the HI 
emission at 5 km s~^ is not exceptionally high, resulting 
in a small ghost correction. After all these tests, we re- 
main confident that the detected variability in the case of 
B1929-hl0 is real. 

5.1.6. PSR B2016+28 (FigureU^ 

In Figure [TOl we compare all four absorption profiles for 
this pulsar. This is the slowest-moving pulsar in our sam- 
ple, so the range of spatial scales covered with our obser- 
vations is only 1-10 AUs. No significant variations were 
found. Our 2-ct sensitivity is A(///o) ^ 0.15 for all base- 
lines. Frail et al. (1994) found A(///o) - 0.15 (from their 
Fig. 3) over periods of 0.6 and 1.7 yr. As first noted in Sta- 
nimirovic et al. (2003), we compared one of our absorption 
profiles for this pulsar with data from Frail et al. and con- 
cluded that slight calibration problems may have affected 
Frail et al.'s data. Two recent pulsar experiments found 



variability at this level but at much larger scales (> 200 
AU; Johnston et al. 2003 and Weisberg et al. 2008). 

5.2. Additional Opacity Fluctuations in Unsmoothed 
Difference Spectra? 

IWeisberg fc Stanimirovid ()2007D presented preliminary 
results of this experiment, and noted that the unsmoothed 
HI absorption spectra exhibited several apparently signifi- 
cant optical depth variations, mostly in the form of single- 
channel fluctuations. The most prominent ones were in the 
direction of B1929-I-10, but there were a few additional 
variable features toward other pulsars. For the current 
study, however, where we have applied hanning smooth- 
ing to ensure full independence between velocity channels, 
only the detections in the direction of B1929-I-10 survive. 
Could the additional variations in the unsmoothed spectra 
be real? 

The minimum expected t emperature for the CNM i n 
thermal equilibrium is 15K ()SpitzeiJ [19781: iHeilesI [2007al ) . 
To achieve such a low temperature, dust grain heating of 
the HI would need to be absent. At this temperature, the 
full-width, half-maximum (FWHM) of an HI feature will 
be 0.8 km s~^ even in the absence of turbulent broad- 
ening; for a less extreme 30K the corresponding width 
is 1.2 km s~^. Using a large sample of CNM clouds, 
iHeiles fc Trolandl (|2003b[ ) deduced a relationship between 
Tk, the temperature inferred from the measured linewidth, 
and the spin temperature Ts : Tk ^ Ts{l + AP-'^^), where 
M is the sonic Mach number. For 196 out of 202 CNM 
components they found that M > 1, implying then that 
Tk/Ts > 1.2. Therefore, even with the minimum expected 
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Fig. 10. — Difference spectra for B2016+28 derived from four different observing epochs. See Fig. [5] for further details. 



turbulence, the HI absorption line FWHM would be ^ 0.9 
km s"^ at Ts - 15 K, or 1.3 k m s~^ at T^ ^ 30 K. So 
while we have rejected most of IWeisberg fc Stanimirovig 
(J2007l )'s variations because unsmoothed data can be mis- 
leading, there are also physical reasons to doubt reality of 
features with FWHM< 1 km s^^. In contrast, our ban- 
ning smoothed resolution appears to match the narrowest 
plausible TSAS features. 

6. DECOMPOSITION INTO GAUSSIAN COMPONENTS 

In this section we derive the spin temperature and the 
CNM column density for gas seen in absorption in the 
direction o f our pulsars. W e emplo y the technique de- 
veloped by iHeiles fc Trolandl ()2003al ). which is based on 
fitting Gaussian components to both HI absorption and 
emission spectra. The purpose of this section is twofold: 
(i) our further data analysis requires temperatures and 
column densities of associated CNM clouds; and (ii) since 
the CNM absorption profiles are often superpositions of 
several components, our original hope was to search for 
temporal changes in both the peak optical depth and ve- 
locity centroids of all individual velocity components. Un- 
fortunately, since the Gaussian decomposition of absorp- 
tio n and emission p r ofiles i s not unique (see discussion 
in iHeiles fc Trolandl (|2003al )). our second goal was not 
achievable. 

To analyze each pulsar, we first selected the epoch with 
the lowest noise level. The HI absorption and emission 
spectra were then decomposed into several Gaussian pro- 
files. We note that three of our six pulsars are at high 
Galactic latitudes and have simple line profiles, easily fit- 
ted with one or two Gaussian functions. Optical depth 
profiles for low latitude pulsars B0540-t-23, B1929-t-10, and 
B2016-I-28, (at h = -3.32, -3.88, and -3.98°, respec- 
tively), are more complex but still with several obvious 



velocity components. 

6.1. Procedure for estimating the spin temperature Tg 

We summarize here very briefiy the basic definitions and 
the procedure, see Heiles fc Troland (2003a) for further 
details. For each pulsar we have the HI emission profile, 
T^(v) (note that this is the emission spectrum that would 
be observed in the absence of the pulsar) and the HI op- 



tical depth spectrum, I{v)/Io = e 



-t{v) 



Both spectra in 



general have multiple velocity components. 

We first fit t{v) with a set of N Gaussian functions using 
a least-squares technique: 

N-l 

r{v) - yZ ro,„e-[("--°-)/^""l' (4) 







where ro,„ is the peak optical depth, vq.u is the central ve- 
locity, and Svn is the 1/e width of component n. N is the 
minimum number of components necessary to make the 
residuals of the fit smaller or comparable to the estimated 
noise level of I{v)/Iq. 

While the optical depth spectrum predominantly re- 
fiects the CNM, both the cold and warm neutral media 
(WNM) contribute to the HI emission spectrum. Conse- 
quently we can write the amplitude of the emission spec- 
trum, Tb(w), as follows: 

Tb(u) = TB.cNAiiv) + Tb,wnm{v). (5) 

We now write each of the two terms of Eq. [S] in detail, 
using the optical depth profile derived from Eq. 2] 

First, Tb,cnm{v), the HI emission originating from N 
CNM components, can be written as: 

Af-l 



Tb,cnm{v) = ^ Ts,n{i 



-T„H)g-E, 



Av) 



(6) 



where T^^n is the spin temperature of cloud n, and the 
subscript m represents one of the M CNM clouds that lie 
in front of cloud n. 
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Next, Ts.wNA'iiv), the HI emission originating from the 
WNM, is represented with a set of K Gaussian functions. 
The complicating factor here is that a certain fraction F of 
the WNM is located in front of the CNM, while a fraction 
(1-F) of the WNM is beyond the CNM with its emissions 
being absorbed by CNM clouds: 

K-l 

Tb,wnm{v) = V[^fc + (l-^fe)e-^(")]xro,fce-[(''-"°''=)/*"''l' 







(7) 

where the subscript k corresponds to each of the WNM 
components and a fraction Fk of the WNM cloud k lies in 
front of all CNM components, while a fraction \ — F^ is 
being absorbed by the CNM clouds. For a given order of 
CNM clouds along the line of sight and a given set of Fk 
values, the T-b{v) profile is simultaneously fitted for Gaus- 
sian parameters of the WNM components and the spin 
temperature of individual CNM clouds. 

For each pulsar, we vary the order of Gaussian functions 
along the line-of-sight (for N CNM components there are 
A^! possible orderings) and perform the Tb(w) fit. We then 
choose the ordering of CNM components that gives the 
smallest residuals in the least-squares fit. Unfortunately, 
the difference in the fit residuals is often not sufficiently 
statistically significant to distinguish between different val- 
ues of Fk- However, Fk has a large effect on the derived 
spin temperatures. Hence we follow the Heiles & Troland 
(2003a) suggestion and estimate the final spin tempera- 
tures by assigning characteristic values of 0, 0.5, or 1 to 
each Fk (among the extreme possible values of and 1), 
and repeating this for all possible combinations of WNM 
clouds. The final spin temperatures are then derived as a 
weighted average over all trials. 

6.2. Spin temperature and HI column density 

Figure [TT] shows results from the Gaussian decompo- 
sition and fitting of HI emission and absorption spec- 
tra for two example pulsars, B1737-I-13 and B1929-I-10. 
This figure has three panels: the top panel, displaying the 
H I emission spectrum, Tb{v), and its decomposition into 
separate CNM [Tb.cnm{v), shown as a dotted line] and 
WNM [Tb,wnm{v), shown as a dashed line] components; 
the middle panel giving the optical depth profile and its 
individual Gaussian-fitted components (shown as dotted 
lines); and the bottom panel showing the difference be- 
tween the observed optical depth profile and its Gaussian 
component fit (dashed lines on this plot show ±lcr noise 
levels). In Table [3] we list the fitted CNM properties. 

With the exception of the B1929-I-10 line of sight, 
all derived spin temperatures listed in Table [3l are in 
agreement with what is t ypically found for CNM clouds 
(iHeiles fc Trolandl 120 03^ Pi ckev et"aIll2003D . For exam- 
ple. iHeiles fc Troland (2003bl) show histograms of WNM 
and CNM temperatures from their HI emission/absorption 
survey of 79 continuum sources. Most CNM clouds have a 
spin temperature in the range 20-70 K, while CNM clouds 
with a temperature in the range 150-200 K appear to be 
rare. However, spin temperatures we derived for CNM 
clouds in the direction of B1929+10 arc high, ^ 150 and 
^ 200 K. This is the first indication that this particular 
line of sight could be sampling atypical ISM conditions. 

In Table S] we provide the HI column density in WNM 
and CNM phases in the direction of each pulsar. The 



first three pulsars have the CNM/total column density of 
about 15-20%, but for B1929-f 10 this ratio is only 7%, 
and for B2016+28 the ratio is about 40%. The very low 
CNM abundance in the direction o f B1929 -I-10 is unusual. 
For comparison, IHeiles fc Trolandl ()2003bD found that all 
except one of their sources with detected CNM have a 
CNM fraction > 10%, with a majority of sightlines hav- 
ing r^ 30%. However, the same study also noted a separate 
class of objects not exhibiting evidence for the CNM along 
their lines of sight. They suggested that at least some of 
these regions may be affected by supershells. The CNM 
fraction in the direction of B1929-t-10 is similar to several 
dir ections with a CNM f r action of o nly 2-4% observed 
b vJStanimirovic fc HeilesI (|2005D and IStanimirovic et all 
(20(r3). Clearly, properties of the CNM in the direction 
of B1929-I-10 stand out as unusual. As we will see in 33 
a significant fraction of this line of sight is located inside 
the Local Bubble. 

7. A STUDY OF THE LINES OF SIGHT 

We have examined the lines of sight of all pulsars in 
our sample. While B1737-M3 and B2016+28 are at a 
significant distance from the Sun and their lines of sight 
probably average over many interstellar clouds, B0823-I-26, 
Bl 133-1-16 and B1929-I-10 are close with significant por- 
tions of their lines of sight probing the very local ISM. We 
detail here some interesting properties of the latter three 
pulsars' lines of sight. 

PSR B1929-I -10 is the closest pul sar in our sample, 
D = 36ltlo pc (jChatteriee et al.ll2004D , and this relatively 
close distance raises the question of whether the observed 
fiuctuations of HI optical depth are somehow related to 
the pulsar's proximity to the solar system. For example, 
we are embedded near the middle of a low-density, high- 
temperature region called the Local Cavity which extends 
1-200 pc before being bounded by a high density wall in 
most directions, so a significant portion of the line of sight 
lies inside it. In what follows, the "Loc al Bubble (LB)" 
consis ts of the Local Cavity plus the wall. iLallement et al.l 
(|2003t ) mapped the dense neutral gas wall bounding the 
Local Cavity via interstellar absorption measurements of 
stars of known distance. The update d map of the neutral 
gas surrounding the Local Cavity bv lWelsh et al.l ()2010a|) 
confirms that the line of sight from B1929-I-10 enters the 
wall of the LB at a distance of ~ 120 pc and then pierces 
the wall for ~ 50 — 6 pc. A lower- r esolut ion interstellar 
reddening map from iToscano et al.l ()1999D also indicates 
that the LB boundary is distended significantly toward 
the pulsar. Furthermore, the B19204-10 line of sight ap- 
pears to be in or near a long extended neutral finger of the 
LB (possibly an interstellar tunnel gouged by a supernova) 
and may therefore lie in the finger or along its wall for some 
distance. The dense finger is most likely associated with 
the Ophiuchus molecular cloud and Ide Geus et al.l (|1990t) 
found that the CO distribution in this region is filamentary 
and has radial velocities in the range 4.2-5.3 km s~^(note 
that our varying HI component is centered at 4.8 ±0.1 
km s~^). 

In addition to the above, we found that two stars within 
3 degrees from B1929-1-10 have been observed in Nal: 
HD178125 is at a distance of 173 pc, while HD180555 is at 
a distance of only 106 pc. Based on the data presented in 
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PSR 



Table 3 
cnm components from gaussian decomposition. 



TO 



vo (LSR) 
(km s"-'^) 



Av 
(km s^"'") 



-^ spin 

(K) 



N(HI)cNM 

(1020 cm-2) 



B0823+26 














0.26 ±0.01 


4.91 ±0.02 


2.3 ±0.1 


62 ±5 


0.72 ±0.02 


B1133+16 














0.16 ±0.02 


-2.10 ±0.03 


1.7±0.1 


27 ±5 


0.15 ±0.01 




0.16 ±0.01 


-3.5 ±0.1 


3.6 ±0.1 


48 ±6 


0.53 ±0.03 


B1737+13 














0.05 ±0.01 


0.7 ±0.6 


4.2 ±0.8 


59 ±6 


0.25 ±0.03 




0.78 ±0.13 


4.19 ±0.06 


2.1 ±0.1 


35 ±7 


1.1 ±0.1 




0.66 ±0.13 


4.59 ±0.05 


1.2 ±0.1 


22 ±6 


0.33 ± 0.05 


B1929+10 














0.017 ±0.002 


-0.7 ±0.3 


5.3 ±0.8 


199 ± 25 


0.35 ±0.04 




0.120 ±0.005 


4.8 ±0.1 


3.0 ±0.1 


148 ±4 


1.00 ±0.04 




0.076 ±0.003 


8.4 ±0.2 


3.9 ±0.3 


206 ±7 


1.15 ±0.06 


B2016+28 














1.12 ±0.03 


3.70 ±0.03 


2.36 ±0.06 


50 ±30 


2.8 ±0.3 




0.67 ±0.04 


5.97 ±0.04 


1.4 ±0.1 


50 ±30 


0.9 ±0.1 




0.81 ±0.05 


9.16 ±0.06 


10.2 ±0.2 


80 ±50 


13.3 ±1.5 




0.57 ±0.05 


9.70 ±0.04 


2.7 ±0.2 


50 ±30 


1.7 ±0.2 




1.30 ±0.03 


13.3 ±0.01 


1.54 ±0.03 


20 ±10 


0.7 ±0.1 



Table 4 
HI Column densities in CNM and WNM toward pulsars 



PSR 


iV(HI)wNM 


iV(HI)cNM 


Ar(HI)TOT 


iV(HI)( 


:;nm/^(HI)tot 




(102" cm-2) 


(102" cm-2) 


(1020 cm-2) 






B0823+26 


4.3 


0.7 


5.0 




0.14 


B1133+16 


3.5 


0.7 


4.2 




0.16 


B1737+13 


6.3 


1.7 


8.0 




0.21 


B1929+10 


34.7 


2.5 


37.2 




0.07 


B2016+28 


28.9 


18.7 


47.7 




0.39 
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iGenova et al.l (|1997[ ). these stars show Nal absorption at 
an LSR velocity of +5.9 and +5.6 km s~^ , respectively (ve- 
locity uncertainty for these measurements is 1-1.5 km s~^). 
As our varying pulsar absorption component is centered at 
+4.8 km s~i (with a FWHM of 3 km s"^), it is most likely 
tracing the same interstellar cloud seen in the direction of 
the two HD stars. This would imply that the absorbing 
cloud is at a distance of < 106 pc, as it is seen against 
HD180555, and therefore embedded in the Local Cavity 
or the LB wall. 

The Local Cavity itself is a region of low scattering; 
while most investigations have concluded that there is 
enhanced scatter ing at its boundary, implying enhanced 
turbulence th ere (iPhillips fc Clegidll992tlBhat et al.l ll998t 
iGupta et all 119991: IChatteriee et al.l |200lD Based on 
modeling of many Nal absorption features, 'Geno va et al.l 
([1997) found that the +5.9 and +5.6 km s^"'^ absorption 
components in the direction of HD178125 and HD180555 
belong to a plane-parallel flow immersed in the Local 
Cavity, the so-called "flow C". They suggested that flow 
C could represent gas condensing out of the LB wall. 
Generally, many diffuse interstellar clouds (so-called "Lo- 
cal Fluff") have been found inside the Local Cavity it- 
self and have been interpreted as signatures of the de- 
velopment of Rayleigh- Taylor instabilities in the wall, 
whereby gas is condensing out of the leading edge of the 
LB wall and flowing toward the Sun (Miiller et al.ll2006t 
iBreitschwerdt fc de Avillej|2006[) . 

To examine lines of sight tow ard B0823+26 and 
B1133+16, we turn to Figure 7 of La llement et al.l (|2003[ ) 
which shows neutral gas in vertical planes at fixed galac- 
tic longitudes. In the case of B0823+26 (b ~ 32°), about 
one-half of the line of sight appears to be inside the Local 
Cavity and traverses the LB wall over only a short dis- 
tance. In the case of B1133+16 {b ~ 69°), the line of sight 
appears to exit the LB through an opening or a chimney 
of low density material rather than a bounding wall or any 
other substantially dense neutral material. 

In summary, since almost one-half of the line of sight 
toward B1929+10 is associated with the LB and the ab- 
sorbing gas is likely to be within 106 pc, we hypothesize 
that our detected TSAS clouds along this line of sight are 
formed or influenced by the LB. Based on the stellar Nal 
absorption lines at similar velocities to those we find in 
HI toward B1929+10, we conclude that our TSAS fluctu- 
ations are likely to be associated with the LB wall. The 
high spin temperature of the absorbing clouds, as shown 
in the previous section, may result from the thermally un- 
stable gas due to wall fragmentation (see also Section 9.4). 
It is then not surprising that we did not detect TSAS in 
the direction of B0823+26 and B1133+16, despite equally 
good sensitivity, because these lines of sight do not tra- 
verse much or any of the wall. 

8. TSAS PROPERTIES INFERRED FROM OBSERVATIONS 

In Sections Uj and |5| we have identified several signifi- 
cant variations in absorption profiles that could be inter- 
preted as TSAS. All of these features are in the direction 
of B1929+10. We summarize their observed properties in 
Table O We measure TSAS properties from EW varia- 
tions, as this is more robust than using difference spectra 
which have a lower S/N. We estimate TSAS properties 



from EW variations with ^ 2.5-cr. 

Table |5| lists the following quantities: the time baseline 
over which the change in EW is observed, the correspond- 
ing transverse spatial scale, the absolute value of AEW, 
the significance of the fluctuation, the HI column density 
of the TSAS (7V(HI)tsas), the inferred HI volume den- 
sity of TSAS (assuming spherical HI blobs whose radius 
is equal I), and the inferred thermal pressure. Some com- 
ments are necessary for interpreting these tables, (i) Our 
observations do not constrain the exact size of TSAS, but 
rather its characteristic scale. In measuring the tranverse 
distance (I) traveled by the pulsar over a time span exhibit- 
ing a variation in HI absorption spectra, we are sampling 
a spatial scale over which we have sampled the optical 
depth function. Therefore, even in the case that signifl- 
cant fluctuations are due to discrete blobs of gas, I does 
not necessarily correspond to a physical size of these blobs. 
This point has been emphasized by Deshpande (2000). 

(ii) The temperature of TSAS, needed to estimate 
iV(HI)TSASj is also not constrained directly by observa- 
tions. We have estimated in Section 6 the temperature of 
the parent CNM clouds with which TSAS is associated. 
In the scenario where TSAS results from discrete, elon - 
gated objects, such as in lHeilesI (|1997^ and lHeilesI (|2007b[ ). 
CNM clouds are constituted from TSAS and an inter- 
TSAS medium. The derived spin temperature in this case 
corresponds to the inter-TSAS medium, while TSAS could 
have a much colder temperature of ~ 15 — 30 K (Heiles 
1997). 

There are several interesting results: 

1. Significant fluctuations (Table [5]) are found on a 
variety of spatial scales I = ^^™ ^ x (6,12,28,46) 
AU (with Z^Hi being a TSAS distance), suggesting 
that the CNM in the direction of B1929+10 has 
structure on spatial scales of 5-50 AUs. Our char- 
acteristic scales for TSAS are similar to results from 
VLBA experiments: Brogan et al. (2005) found a 
typical size of 25 AU and Lazio et al (2009) found 
a typical size of 10 AU. 

2. If we assume that TSAS has the same spin tem- 
perature as its parent CNM cloud {^ 170 K in the 
case of B1929+10), we estimate a HI column den- 
sity of 7V(HI)tsas = (1 - 3) X IQi^ cm'^. This is 
slightly lower than the typical CNM column density 
of 5 X 10^9 cm-2 found for the CNM in the Heiles 
& Troland (2003b) survey. 

3. All TSAS features contribute 10-30 percent to the 
HI column density of their 'parent' CNM clouds. 

4. If we interpret TSAS as being in the form of spher- 
ically symmetric blobs, its observationally inferred 
HI volume density and thermal pressure are > 10'* 
cm""^ and > a few x 10^ cm^'^ K, respectively. 
These values are more than two orders of magni- 
tude larger than what is typically expected for the 
ISM, in agreement with previous TSAS measure- 
ments (see Section 2). 

5. In terms of HI mass, TSAS contains only ^ 10"'^ 
Earth masses. This is many orders of magnitude 
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Fig. 11. — Left: Fitted results for the opacity profile and H I emission for B1737+13 observed at epoch 2001.86. Right: Fitted results for 
the opacity profile and H I emission for B1929+10 observed at epoch 2000.95. (Top) The H I emission profile is shown with the solid line. 
The total H I emission from WNM components is shown with the dot-dot-dashed line. The dashed line shows the emission contribution from 
all CNM components {Tb.cnm(i'))- {Middle) The optical depth profile is shown with the solid line. Fitted Gaussian components are shown 
with dashed lines. {Bottom) Residuals after fitting Gaussian function(s) to the optical depth profile. The overlaid contours show itl-cr levels. 



below what is required for gravitationally confined 
clouds. 

6. In the shaped TSAS model (Heiles et al. 1997) 
TSAS could be much colder with a temperature up 
to a factor of 10 lower than what we have used. This 
would reduce the derived iV(HI)TSAS and ?i(HI) by 
a factor of 10, and P/k by a factor of 100. If an 
elongation factor G ~ 5 — 10 is further introduced, 
the derived thermal pressure can be brought down 
to P/k ~ 2 X 10^ cm~'^ K which is close to the tra- 
ditional hydrostatic pressure (with some turbulent 
contribution). However, this could work only for 
the case of both low TSAS temperature and signifi- 
cant line-of-sight elongation. 

In Table [S] we provide upper limits on TSAS proper- 
ties based on non-detections in the directions of four other 
pulsars. 

9. DISCUSSION 

9.1. The Heiles (1997) discrete cloud model 

In the Heiles (1997) scenario, a typical CNM cloud con- 
sists of an inter-TSAS medium and many discrete TSAS 
features in the form of cylinders or disks, which are homo- 
geneously and isotropically distributed within the cloud. 
An arbitrary line of sight samples many cylinders or disks; 
however only cylinders seen end-on or disks seen edge-on 
have sufficient column densities to produce absorption fea- 
tures recognized as TSAS. Meanwhile, non-end-on cylin- 
ders or non-edge-on disks have very little column density 
and blend to form a smooth background which contributes 
very little to the measured HI absorption profile. The 
model assumed the following basic parameters for TSAS: 
a cloud size perpendicular to the line of sight L± — 30 
AU, Ts = 15 K, iV(HI)TSAS = 0.33 x lO^^ cm^^ and 
K = 3 for the number of end-on cylinders or edge-on 
disks along every line of sight. Two values were consid- 
ered for r, the ratio of the CNM to TSAS column density: 
r = 7V(HI)cnm/A^(HI)tsas = 30 and 100. The model pre- 
dicts the elongation factor G and the volume filling factor 
(j) of TSAS required to ameliorate the over-pressure prob- 



lem, under the assumption that the asymmetric TSAS fea- 
tures are uniformly and isotropically distributed through- 
out CNM clouds. If TSAS is in the form of edge-on disks 
then G ~ 10 and they fill about 3.3% of the CNM volume; 
whereas if TSAS is in the form of end-on cylinders then 
G ~ 3.8 and they ffll about 3.6% of the CNM volume. 
However, the volume filling factor of TSAS is directly pro- 
portional to the number of end-on cylinders or edge-on 
disks H along any line of sight. 

If we assume Ts = 15 K for our TSAS detections, a geo- 
metric factor G = 5 — 10 is required to explain the observed 
TSAS volume density and pressure within the context of 
the Heiles model. This is in the range proposed by Heiles 
(1997). In addition, the TSAS HI column density as well 
as the r ratio are similar to values explored in the model. 
It is therefore reasonable to compare our estimated volume 
filling factor of TSAS with the model predicted values. 

In the case of B1929-f 10, as the pulsar moves through 
the ISM our HI absorption spectra have essentially sam- 
pled a one-dimensional cut across the absorbing CNM 
cloud (about 50 AU long) and detected at least two TSAS 
features. This means that a one-dimensional covering frac- 
tion of TSAS in this direction is at least 40%, if we assume 
a lOAU planc-of-sky TSAS size. We can now estimate the 
fraction of the line of sight occupied by TSAS. We ffi'st es- 
timate the line-of-sight size of the parent CNM cloud. As 
this cloud has N{HI) ~ 10^° cm'^ (from Tabled, and if 
we assume a CNM volume density of 100 cm^"^, this results 
in a length of about 7 x 10^ AU. Assuming a cylindrical 
TSAS with a lOAU plane-of-sky size and a geometrical 
elongation G = 10, results in the volume filling factor of 
< 0.2% in the direction of B1929-t-10. The volume fill- 
ing factor in the case of other pulsars in our experiment 
is zero. We note that our sensitivity for B0823-I-26 and 
B1133-I-16 is equally good as for B1929-I-10, and all three 
pulsars probe interstellar environments with a similar op- 
tical depth (peak of ^ 0.2, see Figures 1 and 2) and a 
simila r line-of-sight length . In a similarly sensitive experi- 
ment, iMi^erleFI^ ()2005[ ) had 18 observing epochs for one 
pulsar, about 150 trials, and no detections. 

The observationally inferred volume filling factor of 
TSAS in the direction of our pulsars (0 or < 0.2%) is lower 
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than the model prediction (~ 3%). More importantly, 
we observe significant inhomogeneity in the galactic dis- 
tribution of TSAS: no detections toward four pulsars, and 
a high detection rate toward B1929-I-10. For the Heiles 
(1997) model to remain viable, our observations require 
inclusion of the possibility of different TSAS properties 
(specifically ij^, G and H) along different lines of sight, 
and/or additional local physical processes which can mod- 
ify the TSAS production rate. For example, a highly elon- 
gated TSAS with L± a. million times smaller than the size 
of its parent CNM cloud would theoretically result in a 
volume filling factor close to zero. While in principle this 
scenario could work, we are still left with the puzzle of 
what physical processes can produce such tiny clouds and 
at the same time greatly vary cloud properties across the 
ISM. 

9.2. The shocked CNM model XHennebelle & Audiil2dd^ } 

Several recent numerical simulations of the small-scale 
structure in the ISM have addressed the TSAS phe- 
nomenon. In particular, iHennebelle fc Audit! (|2007t ) per- 
formed relatively high resolution numerical simulations of 
colliding WNM flows and the formation of CNM clouds. 
A collision of incoming WNM streams creates a thermally 
unstable region of higher density and pressure but lower 
temperature, which further fragments into cold structures. 
The thermally unstable gas has filamentary morphology 
and its fragmentation into cold clouds is promoted and 
controlled by turbulence. While typical CNM clouds are 
formed in these simulations within ~ 1 Myr, transient 
shocked CNM regions, produced by supersonic collisions 
between CNM fragments, are also found, with properties 
similar to those of TSAS: n - lO^"" cm'^, T ~ 30 K, 
P/k ^ 10^ K cm~^, a nd a slightly larger size of ^ 400 
AU. IHennebelle et al.l ()2007D found that 10% of lines of 
sight cross dense CNM gas with n > 10'^ cm^^^, and only 
0.1% of lines of sight cross gas denser than n > 10'* cm~'^, 
suggesting that extremely dense TSAS should rarely be 
detected. 

While this was a 2D simulation of a coherent region 
only 20 x 20 pc large, basic observed properties of the 
CNM/WNM were surprisingly well captured. For ex- 
ample, the CNM column density peak, and the HI col- 
umn density distribution function appear close to observed 
properties based on the Heiles & Troland (2003) survey. 
Considering that the simulation has a resolution of 400 
AU (cruder than the scales generally sampled in pulsar 
and interferometer experiments) and probes regions only 
20 pc across (much smaller than our lines of sight), direct 
comparison with our observations is not straightforward. 
Nevertheless, the predicted properties of the shocked CNM 
clouds, especially cloud abundance, are encouraging. Con- 
sidering that the CNM/TSAS production in converging 
flows is controlled by the level of turbulence in colliding 
flows, local turbulent enhancements may be able to pro- 
vide the necessary tuning and therefore explain the inho- 
mogeneity in the spatial distribution of TSAS observed in 
the direction of B1929+10. 

9.3. The cut-off scale of the CNM turbulence? 

In an attempt to compare various observations, we have 
investigated HI TSAS detections and upper limits from 



all published and recent pulsar and interferometric exper- 
iments as a function of their spatial scale. The results are 
summarized in Figures [12] and 1131 

In Figure [12] we first show the maximum variation 
in HI optical depth (At) as a function of peak optical 
depth along the line of sight (rpcak)- This plot separates 
all sources cleanly into two groups: "low optical depth" 
with Tpcak < 0.5 (B1929+10, B0823+16, B1133-hl6, and 
B0301-M9), and "high optical depth" with Tpcak > 0.5 (3 
interferometric detections, pulsars B1737-I-13 B2016-I-28 
from this study, highly sampled pulsar B03294-54 from 
Minter et al. 2005, and B0736-40, B1451-68, B1557-50 
from Johnston et al. 2003). The only exception in this 
figure is the 2.6-cr detection by Weisberg et al. (2008) in 
the case of B0301-I-19 which probes a low optical depth 
CNM but has At ^0.1. This plot suggests that regions 
with high optical depth (and therefore high CNM column 
density) have higher level of At variations than regions 
with low CNM column density. In other words, this im- 
plies that the the level of t fluctuations tracing TSAS is 
intimately associated with the larger scale structure of the 
parental CNM. 

Figure [13] shows the detected level of optical depth vari- 
ations, or upper limits, as a function of TSAS spatial scale 
/. No obvious correlation between the two quantities is no- 
ticea ble. The Ke ndall's tau correlation test for censored 
data (jlsobe et al.l[l986) gives a probability of 62% that the 
correlation is not present. The low and high Tpcak groups 
of sources nicely separate in this plot around At = 0.03: 
the high- Tpcak sources are above this line, while all low- 
''poak sources except B0301 are below this line. We inves- 
tigated whether a correlation between TSAS scale I and 
At exists for the two groups separately. For both groups 
the Kendall's tau correlation test for censored data gives 
a probability of ^ 50% that the correlation is not present. 

We therefore conclude that the detected level of opti- 
cal depth variations, as well as upper limits, do not show 
evidence for a correlation between I and At. Such a corre- 
lation, with At oc Z("~2)/2^ jg expected in the Deshpande 
(2000) power law model, and would represent the tail-end 
of the turbulent spectrum of the HI optical depth on larger 
scales (0.02 to 4 pc, Deshpande et al. 2000). We overplot 
this predicted level of variations in Figure [13] as a sloping 
line. 

Our compilation of recent interferometric and pulsar 
observations of TSAS does not provide evidence for the 
power-law turbulent spectrum of At on scales < 1 to 1000 
AU. While various sources probe separate lines of sight, we 
note that B0823-I-26, B1133-f 16 and B1929-I-10 all probe 
the local ISM (within 400 pc), with similar peak optical 
depth, and with very similar sensitivity. This essentially 
provides an almost uniformly sampled region on this graph 
from ^ 10 to 200 AU. Similar arguments apply to the 
higher optical depth portion of this diagram. We are left 
with the conclusion that TSAS is a sporadic phenomenon. 

It is interesting to note that the smallest spatial 
scale of TSAS detections in both interferometric and 
B1929-f 10 studies is - 10 AU. Two pulsars (B2016-t-28 
and B0329-I-54) have probed scales < 10 AU, however 
TSAS was not detected. The spatial scales of ~ 10 AU 
for TSAS in the CNM are especially interesting, since in 
the case of hydrodynamic turbulence the cut-off scale of 
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Table 5 
Observed and derived properties of TSAS features with > 2.5 — ct in the direction of B1929+10. 



At I |AEW| Signif. Af(HI)|gAs n(HI)' 



(yr) (AU) (km s'^) (ct) (10 
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'■) (10^ 



') (106 



P/k 



K) 



0.16 6 


0.046 


2.8 


1.4 


16 


27 


0.33 12 


0.083 


3.3 


2.6 


14 


24 


0.75 28 


0.040 


2.5 


1.2 


3 


5 


1.24 46 


0.089 


3.5 


2.8 


4 


7 



"^ We have assumed T^ = 170 K based on Table El 
'' Calculated assuming a spherical geometry. 



Table 6 
Upper limits on observed TSAS properties from non-detections. 



PSR 


I 


|AEW| 


^(HI)^SAS 


n(HI)'' 


P/k 




(AU) 


(km s-i) 


(lOis cm-2) 


(10^ cm-3) 


(10^ cm-3 K) 


B0823+26 


5-50 


0.05 


<6 


<7 


<4 


B1133+16 


20-170 


0.06 


<4 


<2 


< 1 


B1737+13 


20-180 


0.35 


< 19 


<7 


<4 


B2016+28 


1-10 


0.45 


<41 


<270 


< 140 



"^ We have assumed Ts = 50 K. 

'' Calculated assuming a spherical geometry and the largest probed size. 
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Fig. 12. — Variations in HI optical depth (and upper limits) as a function of peak optical depth along the line of sight. Circles (detections) 
and x's with arrows (upper limits) are from this study. Filled squares (detections) and open squares with arrows (non-detections) are from 
Johnston et al. (2003). The horizontal dot-dashed line anchored on either end by arrows represents non-detections over an almost continuous 
sampling of PSR 0329-1-54 by Minter et al. (2005). The triangle at a peak optical depth of ~ 0.3 shows the 2.6-(t detection in the direction 
of PSR B030H-19 by Weisb erg et al.l l(2008|)- Interferometric detecti ons are shown as stars: 3C138 at a peak optical depth of ~ 0.8 (Brogan 
et al. 2005); 3C161, the lower star at a peak optical depth of ~ 1.0 IIGoss et al.ll2008l ): and 3C147, the upper star at a peak optical depth of 
~ 1.0 (Lazio et al. 2009). 
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Fig. 13. — Variations in HI optical depth (and upper limits) as a function of spatial scale. S ymbols are the sa me as in the previous figure. 
The interferometric detections (stars: 3C138 at 25 AU (Brogan ct al. 2005); 3C161 at 500 AU IGoss et al.1[200ll : and 3C147 at 10 AU (Lazio 
et al. 2009). The solid line shows the level of fluctuations in At as predicted by Deshpande (2000). The dot-dashed line represents an almost 
continuous sampling of PSR 0329+54 by Minter et al. (2005) with no detections. 



CNM turbul ence correspon ds to the viscous damping scale 
of ^ 20 AU (|Zweibelll2005l) . However, in the case of mag- 
netohydrodynamic (MHD) turbulence, the turbulent spec- 
trum can cascade below the viscous cut-off scale dow n to 
the scale set by magnetic diffusion (|Cho et al.ll2002l ). In 
fact, magnetic and kinetic spectra can be greatly decou- 
pled, with the magnetic spectrum being much shallower 
than the Kolmogorov spectrum. 

9.4. What's going on with TSAS? 

Our pulsar experiment shows only a few significant 
TSAS detections, all in the direction of one pulsar — 
B1929-I-10. This result strongly suggests inhomogeneity in 
the galactic distribution of TSAS. We have examined two 
proposed origins of TSAS — (i) discrete HI blobs elon- 
gated along the line of sight; and (ii) fluctuations at the 
tail-end of the turbulent CNM spectrum; but our data do 
not provide support for cither of these models. The prin- 
cipal problem is that we do not find these features to be 
homogeneously distributed. Meanwhile, recent numerical 
simulations do generate CNM clouds with a wide range of 
physical properties, including rare examples of very dense 
and over-pressured TSAS-like clouds; but the direct com- 
parison between observations and simulated data is not 
straightforward since simulations probe much smaller ISM 
regions on scales somewhat larger than these experiments. 

The most striking result from our examination of spe- 
cific lines of sight is a possible correlation between the 
B1929-I-10 fluctuations and Nal interstellar clouds inside 
the Local Cavity or at the LB wall. As about one-half of 
B1929-|-10's line of sight is associated with the LB, and a 
significant fraction of this length appears to be traversing 
the LB wall, we suggest that the observed HI absorption 
fluctuations are most likely sampling small-scale structure 
in the LB. Numerous observations have indicated that 
the LB wall is highly inhomogeneous. Many interstellar 



cloudlets found inside the Local Cavity are thought to be 
caused by hydrodynamic instabilities fragmenting the LB 
wall. Through various studies, it has been shown that 
these cloud lets have velocities of up to 25 km s^-^ relative 
to the Sun (jMiiller et al.ll2006l) . This translates to - 5 AU 
yr^^. Coupled with the transverse speed of B1929+10 of 
37 AU yr^^, this suggests that our observed TSAS could 
result from the pulsar's line of sight crossing very local 
cloudlets. 

To investigate the origin of interstellar clouds inside 
the Local Cavity, we look to 2D hydrodynamical sim- 
ulations of the WNM compression by supernova explo- 
sions and subsequent fragmentation into small clouds 
{Koyania fc Inutsuka 2002 ). They showed that within 0.3 
Myr a thermally collapsing layer is formed behind the (su- 
pernova) shock front where cooling dominates heating and 
temperature monotonically decreases. As the gas cools it 
becomes thermally unstable, with a temperature in the 
range of 300-6000 K. Within - 1 Myr this layer starts to 
fragment into many small (< 10^^ pc) cloudlets. Finally, 
these cloudlets can even reach thermal equilibrium at a 
temperature of 20 K and density of 2000 cm"^ . As the ISM 
is frequently compressed by supernova explosions, frag- 
mentation of the shock-compressed medium can provide a 
mechanis m for the formation of many s mall cloudlets. 

In the iKovama fc Inutsukal ()2002t ) simulation, the 
cloudlets become embedded in the hot high-pressure gas as 
the shocked layer fragments, leading to a situation similar 
to what is observed in the LB. The exact te mperature in- 
side the Local Cavity has been uncertain (see I Welsh et al.l 
(|2010bl ) for references), however it is clear that the Cav- 
ity is occupied by a warm/hot ionized medium. The 
cold cloudlets will eventually evaporate in the surround- 
ing warm/hot medium. This process has been studied by 
many autho rs. The mass-loss rate in the case of classic 
evaporation ()McKee fc Cowielll977() is slightly lower when 
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cooling is included (iNagashima et al]|2006[ ). The mass-loss 
rate of lSlavinI (|2007[ ) who considered both flux saturation 
effects and radiative cooling is about 100 times lower than 
the classic evaporation rate. We continue by assuming the 
lowest suggested mass-loss rate, which will have the gen- 
tlest effect on the cold clouds. 

The evaporation timescale depends greatly on the cloud 
size and the temperature of the surrounding medium. For 
a cloud size of 30 AU, Tcvap = 10^ - lO'^ yr, depending 
upon whether the cloud is surrounded by the hot ionized 
medium (HIM) or WNM. For a larger cloud with a size of 
10" AU, Tevap = 10^-^ (HIM) or 10^-^ (WIM) yr. Clearly 
larger clouds survive longer and can even coalesce to build 
larger structures, even if surrounded by a hot medium. In 
about 1 Myr larger clouds could cross a large distance, 
~ 20 pc, and may even escape outside of the LB into the 
ISM through various chimneys and tunnels resulting from 
stellar winds and supernovae, and frequently observed in 
expanding shells. This scenario may be able to explain the 
existence of HI clouds on scales of th ousands of AU or more 
in the ISM, as observed recently by IStanimirovic fc HeilesI 
l|2005[ ). On the contrary, very small clouds will evaporate 
quickly and therefore may not be common in the ISM. 
Their short inferred lifetimes indicate that they would usu- 
ally be found close to their formation sites. 

10. SUMMARY 

We have used the Arecibo telescope to obtain multi- 
epoch HI absorption spectra in the direction of six bright 
pulsars, with the goal of searching for structure on AU 
scales in the CNM. Our target sources were observed previ- 
ously by Frail et al. (1994), who found pervasive tiny-scale 
HI structure. We used an advanced spectrometer capable 
of more accurate measurements in the face of wildly fluc- 
tuating pulsar signals, and carefully analyzed noise statis- 
tics. To search for variability of HI absorption profiles, 
which signifies inhomogeneities on AU scales in the CNM, 
we have searched for time-variability of the absorption 
line equivalent widths and of the absorption spectra them- 
selves. While we had excellent sensitivity, significant vari- 
ations in both EW and absorption spectra have been found 
only in the case of one pulsar: B1929+10. These variations 
imply TSAS on spatial scales I ^ g^^^ x (6, 12, 28, 46) AU 

(where I?hi is the distance to the TSAS feature), with an 
inferred volume density of > 10^ cm""^ and thermal pres- 
sure of > 10^ cm^^ K. 

Our study clearly shows the inhomogeneity in the galac- 
tic distribution of TSAS: the detection rate toward four 
of our pulsars is zero, while the detection rate toward 
B1929-I-10 is high. Such inhomogeneity is hard to explain 
with the model of discrete elongated disks or cylinders 
being a general property of CNM clouds (Hciles 1997); in- 
stead large spatial variations of TSAS physical properties 
or abundances must be invoked. While the Heiles (1997) 



model can help in bringing down the volume density and 
pressure of TSAS, the puzzles regarding its distribution 
and origin still remain. 

Our examination of all pulsar and interferometric TSAS 
detections and upper limits does not show a correlation be- 
tween the level of optical depth fluctuations and the TSAS 
spatial scale, as would be expected if the turbulent spec- 
trum on much larger scales is extrapolated to AU-scales. 
The detections and non-detections probe an almost contin- 
uous range of spatial scales from ~ 0.1 to 1000 AU. There- 
fore, the large number of non-detections of TSAS suggests 
that the CNM clouds on scales 10~^ to 10^ AU are not 
a pervasive property of the ISM. The sporadic TSAS de- 
tections on scales of tens of AU may indicate occasional 
local bursts of turbulent energy dissipation, rather than 
end-points of the turbulent spectrum. 

Another striking result from our study is a possible cor- 
relation between B1929-|-10's TSAS and interstellar clouds 
observed in Nal absorption inside the Local Bubble. We 
find evidence that the TSAS is likely to be within 106 pc 
of the Sun, and is sampling the small-scale structure of the 
LB caused by hydrodynamic instabilities fragmenting the 
LB wall. We propose that the line of sight of B1929-f 10 is 
revealing this recently formed small-scale structure. Simi- 
lar bubbles and their walls are found throughout the Milky 
Way, but the lifetime of a TSAS cloud created from them 
depends strongly on the its size and the temperature of 
the surrounding medium. Larger fragments (size ~ 10'* 
AU) survive longer and can travel large ISM distances, 
becoming a more general ISM property. On the other 
hand, the smallest clouds (size '^ 10 — 100 AU) evapo- 
rate quickly close to their formation site, and are there- 
fore not very commonly observed in the ISM. Additional 
processes, such as stellar mass-loss and collisions of inter- 
stellar clouds/filaments, probably contribute to the CNM 
structure formation on somewhat larger (sub-pc) scales. 
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